Force-induced unfolding of a homopolymer on fractal lattice: 
exact results vs. mean field predictions 
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Abstract 

We study the force-induced unfolding of a homopolymer on the three dimen- 
sional Sierpinski gasket. The polymer is subject to a contact energy between 
nearest neighbour sites not consecutive along the chain and to a stretching 
force. The hierarchical nature of the lattice we consider allows for an exact 
treatment which yields the phase diagram and the critical behaviour. We 
show that for this model mean field predictions are not correct, in particular 
in the exact phase diagram there is not a low temperature reentrance and we 
find that the force induced unfolding transition below the theta temperature 
is second order. 
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The recent development of single molecule techniques have given experimentalists the 
opportunity to grab with suitable handles and mechanically manipulate, by means of optical 
tweezers [[[]] or cantilevers such as atomic force microscopes , proteins, molecular motors or 
DNA molecules. In this way it has become possible to measure or exert on these molecules 
forces in the piconewton range. In particular the effect of a stress on the giant molecule 
titin has been studied in ||-[5| where it was found that a force induced unfolding transition 
takes place between a compact native like state and an extended state. The presence of 
strong hysteresis together with the rather abrupt jumps observed in the force vs. elongation 
curves suggest that the unfolding is a first order phase transition. On the theoretical side, 
many simple models have been proposed to make contact with these experiments (see e.g. 
f6HTT| for the stretching of proteins and homopolymers, |l2|- |l5f| for DNA unzipping and 



|T6Hl8| for RNA unzipping). In particular, in Refs. ||0] the authors study the force-induced 
unfolding of a homopolymer and a heteropolymer in the mean field approximation and find 
that in both cases the critical line separating the globule from the coil is re-entrant at low 
temperatures. This is analogous to what had been found for the phase boundary in the 
phase diagram valid for the DNA unzipping, where it has been proved exactly that in the 



presence of a pulling force mean field is correct [fTJ . 

When self-avoidance is incorporated in the models, exact results are rare and mean field 
treatments become popular. In this work, we study exactly the force-induced unfolding of a 
self-avoiding walk (SAW) on a fractal lattice, the three-dimensional Sierpinski Gasket (SG) 
19| , |20[| . This gasket has a fractal dimension Df = 2, and this renders the system interesting 



from the theoretical point of view because we are below the upper critical dimension for theta 
collapse [pi ]. Mean field theory would predict a reentrance here because the hamiltonian 



compact walks in this lattice have nonzero entropy as compared to the zero entropy of the 
completely stretched coil (same argument as in |14[). However, the exact critical line that 
we find in the SG shows no reentrance. Similarly, by comparing the free energies of the 
globule and that of the coil one would naively expect a first order transition, whereas the 
exact calculation yields a second order transition. Though we cannot conclude from this 



calculation that also on Euclidean lattices the situation will be analogous, we feel that this 
calculation should give a warning that the mean field prediction need not be correct. The 
presence or not of reentrance on Euclidean lattices together with the nature of transition is 
thus an intriguing question which deserves further work. 

We consider a SAW on the three-dimensional SG (3DSG) a hierarchical lattice 

with ramification number 4 and Df = 2 (Fig. 1). We study the combined effect of a force / 
which stretches the polymer along one edge of the SG (Fig. 1) together with a compacting 
self-attractive term, obtaining by assigning a weight exp(/3e) (e > 0) every time two non 
consecutive sites of the SAW are nearest neighbour in the 3DSG. We call (3 = T _1 the 
inverse temperature. To describe the effect of the stretching force, we give each step of the 
walk an orientation and a weight exp(/3/Aa), where / = |/| is the modulus of / and Aa is 
the projection of the oriented step along /. 

The calculation of the partition function and other thermodynamic quantities involves 
the evaluation of twenty-five generating functions (Fig. 2). Twelve generating functions 
involve contributions of a SAW which enters in the n-th order SG at one vertex and goes 
out at one other vertex. Other twelve generating functions arise when an (oriented) SAW 
starts from a vertex, goes out from another one and then re-enters the n-th order SG at a 
later stage. The last generating function is the void. The one-leg generating functions are 
labelled A{, i = 1, ... ,12, and the two- leg contributions B^, i = 1, . . . , 12. One can write 
recursion relations for the generating functions by exactly enumerating all walks in the SG 
on the computer: one has to sum up the contributions of SAWs at the n— th order SG in 
order to generate SAWs in the n + 1-th order SG. 

The initial conditions for the generating functions are: 

Ax = zy' 1 A 2i3 ,4,8 = zy~ 1/2 (1) 

^5,9,10,11 = zy 1/2 A 6J = z A 12 = zy 

#1,2,3,6 = Z 2 W 4 y~ 1 #4,5,8,9 = Z 2 UJ 4 #7,10,11,12 = Z 2 W^y 

where z is the SAW step fugacity, y = exp (/3f), and w = exp (/3e), is the weight responsible 
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for the theta collapse in absence of force. 

Note that, as is apparent from the above equations, we have adopted the convention 
that interactions are restricted to sites within the first order SG and moreover that after 
one SAW has touched one vertex in the SG it is obliged to exit from the SG at that order 
[ 2D| , |22|pB| . It has been proved that this approach is equivalent to the more general approach 



as regards universality of the phase transition |23| . 



When there is no pulling force the recursion relations simplify because: 

A = A Bi = B Vz = l,...,12, (2) 

so one gets |20| : 



A' = A 2 + 2A 3 + 2A 4 + AA 3 B + 6A 2 B 2 , (3) 
B' = A A + AA 3 B + 22B 4 . 

At non-zero /, the equations retain the same structure but every B and every A has 
to be labelled by the appropriate number as found in the exact enumeration. To solve the 
model, we proceed as follows. First, we must find the fixed points of the recursion relations: 
in general, for every phase or critical point (line) there is a fixed point. Then, for every / 
and T fixed, we have to find the critical step fugacity, z c ((3e,/3f) such that the flux defined 
by the twenty-five recursion relations and initial conditions will bring the system in the 
fixed point corresponding to those values of / and T. This critical step fugacity z c (/3e,(3f) 
allows us to have a numerical expression for all the quantities we are interested in, namely, 
the free energy F(f3e,/3f) = 4 log (z c (/3e, /3f)), the average elongation along the direction 
of the pulling force, (x)(/3e, (3f) and the average number of contacts, (n)(/3e, /?/), found by 
calculating the appropriate derivative of F(f3e,f3f). 

From the linearized recursion equations , it is possible to get the critical exponents of 
the phase transition. Every eigenvalue A (there are at most as many such eigenvalues as 
there are recursion relations) , of the linearized flux at given / and T defines two critical 
exponents, Y and u, through: 



A = 2 Y = 2i (4) 



where v = y is defined in terms of the average squared elongation (x 2 ) as (x 2 ) ~ iV 2i/ [21 
for large values of the number of steps in the SAW, N. 

Performing this analysis for our model, we find six different fixed points. At zero force, 
we re-obtain the fixed points given in f2"0|| . In particular, for T < Tg = lo 3 3 , , the theta 
collapse temperature in the absence of force, there is a zero force compact fixed point for 
the recursion relations: 

At = Bi = 22- 1/3 ee B* Vi = 1, . . • , 12. (5) 

For T > Tg, i.e. in the swollen phase, the fixed point that is approached is: 

Ai ~ 0.4294. . . 5,- 0.04998... \/i = 1, . . . , 12. (6) 

Just at criticality at the theta temperature, T = Tg, the fixed point is: 

Ai = 1/3 Bi = 1/3 Vi = 1, . . . , 12. (7) 

The corresponding critical indeces v (corresponding to the largest eigenvelues, Ai, in every 
one of the three regimes) are 1/2, 0.5294... and 0.7294... respectively in the collapsed 
phase, at the theta point and in the swollen phase. 

At T < Tg, there exists a critical line / = / C (T) separating a compact from an open 
phase. When / < / C (T) the recursion relations display another compact fixed point: 

Ai = \/i = 1, . . . , 12 (8) 

-84,5,8,9 = 22 5i, 2,3,6 = -87,10,11,12 = +00 

The divergence arises from the pulling force and one can convince that the two leg generating 
functions are diverging (vanishing) for L — > 00 {L is the system size, and one has L = 2 n_1 
at the n-th level of iteration in the construction of the 3DSG) as y L {y~ L ). Consequently, 
we can eliminate the divergences if e.g. we multiply the diverging functions for Bi, and 
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the vanishing ones for B7, and then take the square root of the result: if we do this all the 
two leg generating functions converge to the fixed point 22 -1 / 3 as in the zero force compact 
phase. When the SAW is in the open phase (/ > f c (T)), the fixed point is: 

A = Vi = l,...,ll; A 12 = l (9) 

Bi = Q Vi = 1, . . . , 12 

Finally, when the force is exacly tuned at the critical line, the fixed point which is approached 
is: 

A i = Vi = l,...,ll; A 13 = rrr ~ 0.5668... = A* (10) 

' ' ij 1 + 6*22- 2 / 3 v ; 

-84,5,8,9 = 22 -1 / 3 -Bi, 2,3,6 = -87,10,11,12 = +00 

again the square root of the products of one diverging two leg function times a vanishing 
two leg function is 22 -1 / 3 . We will see later that it is possible to justify the value found 
numerically also for the nonzero force system. The exponent v is equal to 1/2 for / < f c (T) 
and is 1 in the open phase. 

In practice, by taking advantage of the knowledge of the fixed points, one can devise a 
convenient way to calculate numerically the critical force as a function of temperature. Once 
the critical step fugacity z c {(3e,(3f = 0), of the compact phase at a fixed temperature T, is 
known one can further fix z to this value, and then tune / in order to pass from the compact 
fixed point to the critical line fixed point (within the precision allowed by the computer). 
In order to have a sufficient accuracy in these calculations, it proved necessary to use a 
computer with quadruple precision. 

The phase diagram obtained numerically is shown in Fig. 3. Two remarks are in order. 
First, we notice that there is no re-entrance in the critical line, the slope at T = being 
approximately —0.05. This is at variance with the prediction doable on the basis of mean 
field-like treatments (as those done e.g. in Refs. see also below). Second, the behaviour 
of the critical line near Tg is: 

/ c (T)~(T e -T) a=0 - 87±0 - 01 , (11) 
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in agreement with the prediction a = ^ = 0.868, where v$ is the critical exponent of the end- 
to-end distance at the theta temperature and <fi is the theta transition crossover exponent 
PH (see below for an argument, leading to the behaviour in Eq. |TT] ) . 

We report in Fig. 4 the plot of (x) vs. / at T = 0.35. This figure supports the hypotesis 
that the unfolding transition for the homopolymeric SAW on the 3DSG is second order. 
This is in agreement with the argument based on the recursion relations given below. Our 
belief supported by the exact numerics is that the transition is second order for any nonzero 
T and is first order only at T = 0. 

We now present some arguments to interpret our results (based on the renormalization 
group (RG) flux). Even though to write explicitly the equations and to cope numerically 
with them it was necessary to put the force in step by step, it is equivalent to evolve 
the / = generating functions (Eq. |3[) and then put in the dependence on / at every 
iteration by multiplying the generating functions times suitable powers of y. If we do this, 
we can exploit the symmetries of the problem which make so that before multiplying the 
generating functions times y dependent term, all one leg and two leg diagrams are separately 
equivalent. By noting the structure of the nonzero force fixed points, one can argue that 
the one leg generating functions must vanish in the compact phase like A ~ Cy~ L , where 
y c = exp (/3f c (T)) and C is a constant, for now undetermined. By matching the exponentials 
in Eq.^, we obtain that the two equations, when T <^Tq (T — > 0), can be approximated by: 

A' = A 2 + QA 2 B 2 , B' = 22B 4 . (12) 

The flux corresponding to these RG equations is shown in Fig. 5. The non-trivial fixed point, 
(A*, B*) with both A* and B* non-zero, is obtained with B* = 22~ 1 ^ and A* = - * ,- 2 . 

v ' l+b(B *) 

Thus A ~ Cy~ L as L — > oo, with C = A*. When T — > 0, we notice that the fixed points can 
be approached at once, i.e. since the initial condition, also in the original (not simplified) 
recursion relations, because the terms that should evolve to have an initial value which 
is zero for T — > 0. As a consequence, as T — > 0, the phase boundary of Fig. 3 is found 
by matching the initial conditions with the fixed point values for A& and any B. Thus the 
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critical line is found by solving this system: 

z c = exp(-2[3e)(B*) 1/2 , (13) 
z c = exp (-/3f)A*, 

where B* and A* are the nonzero fixed points of Eq. [12|, corresponding to the critical force 
(such that y L A ~ 1 for large L, see Eq.|10D . The critical line for T — > is: 

f c (T) ~ 2e + T log (^p) ~ 2e - 0.0525 . . . T, (14) 

so that f c (T) starts with negative slope as found numerically. For / > f c (T), the fixed point 
to be reached in the (A, B) plane of Fig. 5 is (1, 0). We observe that for a small deviation 
from A* the RG flux takes the SAW to the open fixed point by moving along the line B(A), 
which we can find in the neighbourhood of (A*, B*) by requiring that it be a fixed line under 
the flux defined by the recursion equations. We require that B(A) ~ B* + c(5A) a in the 
neighbourhood of A*. We thus need to solve the system: 



A' = [A* + (5A)] 2 1 + 6 [B* + c(5A) a ] 2 (15) 



B{A') = 22 [B* + c{5A) a ] 4 

in the unknown quantities c and a. We obtain c = — 1/12A* 3 B* ~ —1.2826 . . . and a = 2 so 
that the fixed line B(A) approaches smoothly the fixed point at critical force. By inserting 
these values in the expression B(A) ~ B* + c(5A) a together with the ansatz: 

z c ((3f c , (3e) - z c (Pf, (3e) ~ (/ - f c )\ (16) 

one gets 7 = 2 which implies (x) ~ (f—fc) f° r / ~ fc and the transition is second order. This 
argument is strictly valid for low T. However it is unlikely that the order of the transition 
could change along the phase boundary and this is also confirmed by our numerics. The 
situation at the point T = is somewhat special: the entropy vanishes and balancing the 
energetic terms gives a first order transition. 

One should notice the importance of the term 6A 2 B 2 . Let us think of a mathematical 
simplified model in which the relevant equations are: 
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A' = A 2 , B' = 22B\ (17) 

in which the mixed term (physically due to stretched walks that still make an extensive 

number of contacts) is suppressed. This corresponds to balancing the free energy of a 

stretched coil with that of a compact globule. In this simplified treatment the re-entrance 

is present, as can be expected from naife estimates of the ground state entropies of the 

stretched and compact state, and the transition is first order. 

Finally, we argue that the exponent a in Eq.|y] is given by a = ^ as anticipated above. 

Near Tg at zero force the free energy behaves as F(T) ~ (Tg — T) 2 ~ a ~^ . For / ~ 0, on 

i 

the other hand, at T = Tg, the free energy behaves as F(f) ~ /"<?. Consequently, we get 
a = ^ in Eq. [0]. This is in agreement with the result found in d = 3 in Ref. || where the 
exponents both take their mean field values (1/2), and also with the mean field analysis in 
Ref. H|7| in which = 1 and consequently a = v (1/2 in the ideal case treated in 0|7|]). 

In conclusion, we presented an exact calculation of the phase diagram of a SAW in 
presence of a compacting contact energy and a stretching force. We deem it is interesting 
because it can be analyzed exactly. A mean field like treatment gives a reentrant boundary 
and a first order transition. Both these predictions are not confirmed by the exact treat- 
ment, which gives a critical line with a zero temperature negative slope and a second order 
transition. This has been explained in a simple way by analyzing a simplified version of the 
recursion relations, analogous to the real space RG equations. Whether or not the critical 
line in the hypercubic lattice shows reentrance is therefore not yet clear and appears to 
be an intriguing question. The 3DSG has Df = 2 and so the most natural comparison is 
with the Monte-Carlo simulations performed in Ref. ||, which indeed gives a second order 
transition, even though we cannot be sure that the order of the transition is the same in the 
two-dimensional real and fractal lattices. The behaviour of the critical force near the theta 
point has been found and a general argument, valid also for hypercubic lattices, has been 
given in agreement with the result found here and also with the result found in d = 3 in [|8] . 

This work was supported by cofin2001. 
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FIGURES 




FIG. 1. In this figure we show the 3DSG at the second stage of iteration construction. Four 
tetrahedra (dashed lines) at the first order of iteration are put together in order to form the 
tetrahedron which constitutes the 3DSG at the second iteration (bold lines). The stretching force 
/ acts along one of the edges of the 3DSG and is also shown in the figure. 
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FIG. 2. In this figure we show the generating functions at the first order together with their 
names given in the text. 
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FIG. 3. Plot of the phase diagram on the 3DSG found numerically. In this figure we have taken 
e = 1/2 to make the calculations. 
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FIG. 4. Plot of the average elongation scaled by N in the thermodynamic limit as a function 
of f for T = 0.35. The critical force found numerically by imposing that the fixed point reached 
after iteration of the recursion is that of Eq. |l^ is approximately 0.969 .... We have taken e = 1/2 
in these calculations. 
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FIG. 5. Plot of the flux lines resulting from the simplified recursion relations corresponding to 
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